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I.  lON/ION  COLLISIONS  AND  MAXWELLUN  EQUILIBRIUM 
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Consider  two-body  collisions  among  monoenergetic  particles  of  density  n^  (1/cm  ), 
charge  Z,  at  energy  £  (eV).  The  collision  frequency  f.  is  taken  as  the  inverse  of  the 
two-body  energy  exchange  (90*  deflection)  or  slowing-<iown  time,  t^  (Delcroixi),  given  by 

=  l/t,  =  1.35E-7  [n_^Z<LN(r)l/[A'’'(E^)*-^  (1) 

Here  A  is  the  ion  mass  number  and  T  is  the  Coulomb  impact  parameter.  The  two— body 
collision  rate  density  is  just  fusion  density  is 

Of  =  (biP(n/)(ajV^)  (2) 

where  b.^  is  0.25  (0.5)  for  differing  (identical)  ions,  the  ratio  of  two-body  to  fusion 
collisions  is  =  q^/q^ 

qj,  =  1.08E11  (Z<LN(r)l/(b.,a^E/]  (3) 

for  collision  energy  (CM  frame)  E^  in  eV  and  a^in  barns.  Here  the  particle  collision  speed 
has  been  taken  to  be  v  =  [2£  /rnkf  ’^  for  monoenergetic  particles  (all  fusion  reactions 
have  been  assumed  to  be  at  90*  incidence)  and  is  to  be  evaluated  at  v^. 

If  ((7fV^)  is  integrated  over  the  actual  collisional  energy  distribution  from  the  initial 
radially-monoenergetic  ion  energy  distribution,  fusions  are  more  probable  and  the 
coefficient  in  eq.  (3)  will  become  considerably  smaller  than  given  above.  For  example,  for 

t 

DT,  if  E  is  taken  as  10  keV,  and  the  initial  collision  energy  distribution 
c 

E^(^)  =  E^(l-cos^)  is  used,  then  the  integrated  average  value  of  <av>  is  about  3.4  times 
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larger  than  for  the  monoenergetic  calculation.  Thus,  initially  -  at  this  mean  energy  -  the 
ratio  q2f  will  0.29  times  that  calculated  from  eq.  (3);  at  =  20  keV  the  factor  is  about 
0.62. 


The  mean  time  for  "approximate  Maxwellianization"*  (t^)  can  be  taken  as  the 
collisionaJ  slowing-down  time  (tj,  from  eq.  (1),  at  the  core  high  density  conditions, 
divided  by  the  fractional  time  (g.)  that  any  given  ion  spends  in  this  region,  thus 

shown  by  Rosenbluth  et  al*  the  time  (tj^)  to  Maxwellianize 
particles  by  upscatter  to  velocity  Vj^  in  the  high  energy  "tail"  of  the  distribution  is  greater 
than  this  (because  of  cross-section  decrease  with  energy)  by  a  factor  of  about 
(l/2)(vj^/v^)  ,  thus  the  time  to  Maxwellianize  particles  to  energy  in  the  "tail"  is 

It  is  evident  that  an  increase  from  E^  to  6E^,  say,  requires  about  7.2t^.  If  E^  is 
initially  at  10  keV,  for  example,  then  the  tail  of  the  distribution  will  be  filled  at  60  keV  in 
this  time.  Now,  at  this  higher  energy  the  fusion  cross-section  for  the  fast  particles  (of  DT) 
with  bulk  ions  is  about  20  times  larger  than  for  monoenergetic  10  keV  ions  alone,  thus  their 
fusion  reaction  rate  will  be  higher  by  about  20(60/2)®’®  «  110  than  that  of  the  (assumed) 
initial  10  keV  population.  These  two  effects  plus  the  increase  in  ion  speed  with  (E^)®'®  are 
responsible  for  the  well-known  fact  that  most  fusion  reactions  in  thermalized  DT  at  ca.  10 

keV  come  from  regions  at  50-60  keV  in  the  Maxwellian  tail.  _ I 

Acrnmsim  Fo» 

VfiC  tAB  j 

*  At  t  =  the  mean  energy  dUtribution  will  have  relaxed  lo  that  the  low—  and  md— range  distribution  rlWtian/ 

is  "approximately  Maxwellian".  •liability 
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For  monoenergetic  ions  at  10  keV,  is  about  q2f  »  3.1E4,  while  at  60  keV, 
q^f  B  210,  corrected  for  all— angle  collisions.  With  such  coUisionality  the  bulk  collisions  in 
the  central  core  (r  <  r^)  will  tend  to  distort  the  core  in  energy  distribution  from  its  initial 
form.  However,  it  is  important  to  note  that  the  two-body  collisions  in  eqs.  (1)  and  (3)  do 
not  involve  the  same  ions  from  one  moment  to  the  next.  This  is  because  each  ion  resides  in 
the  core  only  a  small  fraction  of  its  transit  time  across  the  system,  thus  the  rate  of 
collisional  energy  scattering  is  much  less  than  would  be  adduced  from  the  calculation  of  q2j, 
above. 


The  time  any  given  ion  spends  in  the  core  is  just  6t^  =  while  the  time 
required  for  a  system  transit  is  6tj^  =  2fdr/v(r)  integrated  from  0  to  R.  For  an  assumed 
parabolic  potential  well  shape  the  ion  speed  is  given  crudely  by  v(r)  =  v^[l-<r>^]°‘^, 
where  <r>  =  (r/R).  With  this  the  transit  time  is  estimated  as  6tjj^ «  i^R/v^,  so  that  the 
fraction  of  time  spent  in  the  core  is  given  by  g.  «  tJrR  =  <r^>/'K. 

Thus  the  number  of  upscattering  two— body  collisions  per  ion  per  fusion  reaction 
becomes  q2j{^^  =  q2j<r^>/‘K,  so  that  the  previous  examples  give  only  qj^^^  »  101  for 
=  1E4  eV  and  ~  Maxwellianization  requires  at  least 

t^  for  ions  in  the  bulk  and  (t^/2)(Ej^/E^)^’®  for  the  tail  ions,  it  is  evident  that  full 
Maxwellianization  can  not  occur  for  well  depths  of  25-30  keV  or  greater,  and  that  the  high 
energy  tail  distribution  will  never  be  filled  for  such  mean  core  energies. 

Core  mean  energies  below  this  collisional  energy  range  can  result  in  Maixwellian 

distributions  well  out  into  the  tail.  For  Maxwellian  distributions  in  the  core,  the  qj^  ratio 

can  be  written  for  DT  with  a  Maxwell-averaged  <av>  at  core  mean  energy  E  as 

.  m  m 
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W  =  lI'N(r)l/(«n->(EJ*  '| 


(4) 


f 


Taking  the  Coulomb  factor  to  be  LN(r)  =  10  as  typical,  this  yields  «  200  at 

ts  5£4  eV  as  the  average  number  of  two— body  energy  exchange  collisions  per  average 
fusion  event.  However  here,  as  before,  ions  in  the  tail  can  have  much  higher  fusion  reaction 
rates,  thus  will  tend  to  be  reacted  before  reaching  full  Maxwellian  density  equilibrium  for 
<<  E-at-peak-cross-section. 

Now,  ions  with  energies  Ej^  >  E^  in  the  Maxwellian  tail  will  circulate  to  radial 
portions  further  out  than  the  ion  injection  radius,  <  R.  As  the  Maxwellian  tail  is  built 
up  by  continuous  collisions  within  the  small  central  core,  these  fast  ions  will  climb  further 
and  further  out  until  they  are  eventually  lost  to  the  system.  If  the  well  is  three  times 
deeper  than  the  potential  at  the  injection  radius,  fast  ions  must  have  resided  in  the  core 
collision  system  for  (3^'®)/2  »  2.6  times  the  two-body  collision  time,  t^.  If  the  well  is  six 
times  deeper,  the  fast  ion  residence  time  would  be  ca.  7.2t^.  In  this  case  the  fast  ions  at 
Ej^  >  6E^  will  all  escape  from  the  system.  The  energy  they  carry  at  escape  is  small, 
however,  in  that  their  fast  tail  energy  is  lost  to  the  confining  well  before  their  departure 
bom  the  system. 

Ion  collisional  makeup  will  be  required  to  supply  the  thus-truncated  tail,  but  at  a 
rate  of  only  about  2(E^/Ej^)^’®  «  0.13  times  that  of  the  basic  two-body  collis.on  process. 
This  energy  exchange  rate  is  just  qj  =  nfj,  from  eq.  (1),  and  is  q^  «  1.0E21  coll/sec/cm^  in 
the  core  with  n^  =  lE17/cm*  at  =  5E4  eV.  If  the  core  radius  is  r^  =  1  cm,  the  total 
collision  rate  =  q2(core  volume)  is  Q  »  4.3E21  coll/sec.  However,  any  given  ion  is  only 
in  the  high  density  core  region  for  a  fraction  g.  of  its  lifetime  within  the  system,  thus  the 
rate  of  upscattering  of  each  ion  in  the  system  will  be  less  by  this  factor  than  would  be 
estimated  from  the  expression  for  loss  to  a  sink  at  energy  E^  above  core  energy  E^  (or  E^). 
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Thus  the  actual  loss  rate  for  this  example  would  be 


(5) 

The  ion  in-coie  time  fraction  is  (as  previously  shown)  about  g.  =  <r^>/ir.  With 
this,  for  a  convergence  ratio  of  <r^>  =  lE-2,  the  previous  example  gives  a  DT  mass  loss 
rate  of  (dM/dt)}^^^  »  2.74E-4  gm/sec  to  >  6E^.  This  is  a  volumetric  loss  rate  of  about 
1.5E-3  liter  atm/sec,  adding  a  leakage  gas  load  to  the  vacuum  pumping  system.  Note  that 
the  energy  of  these  escaping  ions  outside  the  system  is  not  6E^  but  is  very  small;  most  of 
their  energy  has  been  given  back  to  the  well;  this  loss  should  not  create  any  significant 
problem  of  ion  sputtering  external  to  the  system.  If  E^  =  2E4  eV,  the  loss  rate  is  four 
times  greater  for  E.  >  6  E  and  eleven  times  greater  for  E^  >  3  E  .  Loss  rates  for  this 

Li  m  L  in 

latter  case  would  be  1.87  mg/sec  or  1.68E-2  liter  atm/sec.  If  pumping  capacity  is  3E4  liter 
atm/sec  this  would  limit  the  system  base  pressure  to  about  5.7E-7  atm. 

The  system  net  energy  gain  (balance)  will  be  little  affected,  since  most  energy  loss  is 
due  to  electron  escape  through  the  cusps.  Note  that  this  gain  IS  reduced  by  the  factor  of 
increased  electron  injection  energy  required  to  create  a  deeper  well  to  hold  upscattered 
ions,  vs.  the  minimum  well  depth  needed  for  core  ion  energy  alone.  But  deeper-than- 
core-energy  wells  are  required  anyway,  to  allow  ion  injection  deep  enough  within  the 
magnetic  field  surface  to  avoid  excessive  transverse  ion  momentum  generation  in  core  ion 
recirculating  flow. 
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n.  ELECTRON  THERMALIZATION 


Electron  collisional  thennalization  will  also  occur,  but  principally  in  regions  of  large 
density;  and  in  these  regions  the  electron  energy  is  low.  In  the  core  the  relevant  collision 
type  is  that  of  ion  collisional  heating  of  electrons.  The  relaxation  time  for  equilibration  in 
such  collisions  is  (m./m^)  longer  than  that  for  electron  self-coUisional  relaxation.  Thus, 
since  ion/ion  two-body  collisional  relaxation  is  slower  by  the  square  root  of  this  factor,  ion 
heating  of  electrons  will  be  slower  than  ion  two-body  collisions  by  this  same  factor, 
(ni./in^)°^ 


There  are  two  points  of  interest  here.  First  is  the  electron  thermal— ization  time 
compared  to  the  electron  lifetime  in  the  high  density  core,  as  limited  by  transits  before 
cusp  escape.  Spitzer^  gives  this  time  as 


f4.0E6  A  A.] 
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For  a  core  radius  of  =  1  cm,  an  assumed  electron  energy  of  e.g.  1  keV  in  the  core, 
so  that  »  1.8E9  cm/sec,  and  an  electron  current  recirculation  factor  of  =  1E5,  the 
core  lifetime  is  t^^  »  I.IE— 4  sec.  At  this  electron  energy,  with  ion  energy  of  E^  =  5E4  eV 

A 

and  LN(r)  =  20,  the  ion/dectron  equilibration  time  for  n.  =  lE17/cm  ,  for  example, 
would  be  t^^  »  6.7E-6  sec  from  eq.  (6),  above.  Thus,  1  keV  dectrons  will  be  wdl 
thermalized  after  only  a  few  thousand  passes  through  the  core.  The  assumed  (low)  electron 
energy  therefore  is  not  valid. 

At  the  other  extreme,  assume  that  the  dectrons  have  energy  of '5E4  eV  (as  though 
already  in  equilibrium  with  core  ions).  Then  the  equilibration  time  is  t^^  s  2.1E-3  sec. 
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while  the  lifetime  in  the  core  is  only  t^^  s  1.6£^  sec;  about  130  times  less  than  this. 

Thus  elections  in  this  example  will  not  remain  cold,  nor  will  they  reach  ion 
temperatures/energies.  Rather,  the  electron  energy  will  tend  to  reach  an  equilibrium 
value  when  t^^  =  t^^,  or  when  the  core  lifetime  2Gjr^/v^  s  1.7E8(E^)^  ®/(n.Z^LN(r)). 

For  V  s  (2E  /m  this  gives  the  equilibrium  core  electron  energy  as 

EJ^^  «  1.45E-8Z[n.r^G.LN(r)]°  ®  (eV).  If  n.  =  lEl7/cm^  r^  =  1  cm,  LN(r)  =  20,  then 

E^{^^  =  6.7E3  eV.  For  higher  ion  density,  e.g.  n.  =  lE18/cm^,  this  is  still  only  «  20  keV. 

However  if  G^  ->  1E6  the  calculated  equilibrium  electron  temperature  for  this  last  case 

becomes  an  unphysical  63  keV,  and  electron  losses  would  become  greatly  enhanced  by  core 

thermalization. 

For  more  realistic,  lesser  (initial  experimental)  conditions,  electron  thermalization  is 
not  significant.  For  example,  for  G^  =  1E3,  n.  =  lE12/cm^,  and  =  1  cm,  this  simple 
model  gives  E^j^^  »  2.1  eV.  Since  this  is  very  much  less  than  (m^/m.)E.  for  all 
E.  >  2E4  eV,  the  ion  collision  rate  with  electrons  would  dominate,  and  the  equilibration 
electron  temperature  due  to  ion/electron  collisions  is  expected  to  be  only 
E^leq  »  (m^/m.)E..  This  does  not  take  into  account  electron  acceleration  by  a  central 
virtual  anode,  which  can  raise  the  electron  energy.  In  general  it  is  expected  that  core 
collisional  equilibration  of  electrons  will  be  insignificant  at  planned  laboratory 
experimental  conditions  and  will  become  noticeable  but  not  dominant  as  reactor  conditions 
are  approached. 

Electron  collisions  will  also  occur  outside  the  core  region;  here  the  electron  density 
will  tend  to  follow  the  ion  density.  Taking  a  parabolic  distribution  as  illustrative, 

“  typical,  where  n^^  is  electron  surface  density  (n^^  =  G^  and  the 

n  ^ 

electron  energy  variation  also  as  parabolic,  E  (r)  s  £  (r/R)  ,  yields  the  functional 
dependence  of  the  equilibration  time  as  t^  e  K(E^^)^'®<r^>®(r/r^)®/n^^,  where  K  is  a 
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factor  contaixiing  design  parameters  and  constants  of  the  system.  From  this  it  is  evident 
that  electron  equilibration  will  proceed  much  more  slowly  as  r  ->  R  than  as  r  ->  r^,  and 
hence  can  be  neglected  at  least  to  first  order. 

in.  lON/ION  COLLISIONAL  DISPERSION 

Returning  to  ion  collisions,  it  is  of  interest  to  estimate  spreading  of 
initially— radially— monoenergetic  energy  distributions  and  of  the  ion  spatial  convergence 
radius  due  to  ion/ion  self-collisions  in  the  core.  Lovberg*  has  shown  that  these  effects  can 
be  kept  non-dominant  in  system  operation. 

The  first  complete  analysis  of  coUisional  equilibration  of  particles  with  initial 
uniform  and  equal  energy  was  carried  out  by  Maxwell^  in  1859,  who  showed  that  the 
stable  state  distribution  of  energy  (or  velodty)  was  (what  we  now  call)  the  Maxwell  (or 
Maxwdl— Boltzmann)  distribution 

fj^  =  (ra/2irkT)^'®EXP(-mvV2kT)  (7) 

In  the  central  core,  ion/ion  collisions  will  produce  this  distribution,  albeit  with  a 
truncated  tail  (for  some  conditions,  as  explained  above)  with  the  consequences  described 
previously.  This  distribution  is  achieved  by  the  combined  effect  of  a  very  large  number  of 
small-angle  collisions  of  each  particle  with  the  Coulomb  fields  of  the  surrounding  particles 
that  can  affect  its  motion.  These  are  limited  to  those  contained  within  a  sphere  around 
each  particle  with  radius  equal  to  the  Debye  shielding  radius 

Ljj  =  (kT/W)®-®  (8) 
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For  ion  core  energies  >  1E3  eV,  the  Debye  sphere  will  contain  a  very  large 
number  of  ions  for  any  density  n  <  lE19/cm^,  thus  most  collisions-deflections  will  be  due 
to  distant  small— angle  encounters.  Any  one  collision  will  yield  only  a  small  deflection  and 
momentum  transfer  to  the  deflected  particle.  The  next  collision  will  be  uncorrelated  in 
direction,  and  successive  collisions  equally  so.  The  net  result  is  that  the  spread  of  particle 
motibn  and  momentum  (and  thus  of  energy)  will  proceed  as  a  collisional  random— walk 
process  until  momentum  equivalent  to  the  particle  initial  momentum  is  exchanged  with  the 
plasma  system. 

At  this  point  the  particle  will  have  been  deflected  approximately  ir/2  (90* )  from  its 
initial  course  and  can  be  viewed  as  then  having  undergone  an  "energy  exchange  collision" 
of  the  sort  used  in  the  preceding  sections  herewith.  While  the  time  to  reach  the  state  of 
this  "collision"  and  the  energy  exchange  associated  with  it  are  both  useful  in  zero  order 
analyses  (as  above)  they  do  not  disclose  details  of  the  random  walk  process  that  are 
important  in  assessing  first  order  effects  in  core  collisionality. 

To  illustrate  these  it  is  useful  to  consider  the  random  walk  problem.  Following 
Jeans*,  the  random  collisional  one-dimensional  motion  of  a  particle  traversing  a  mean  free 
path  £  between  collisions  is  describable  by  the  probability  P  that  it  will  have  moved  a 
distance  s£  in  N  collisions,  by 

P  =  N!/p!q!2^  (9) 

where  p+q  =  N  and  p-q  =  s.  Noting  that  N  is  very  large  in  problems  of  interest  and  that 
8  is  small  relative  to  N  allows  this  factorial  representation  to  be  reduced  by  Stirling’s 
approximation  (for  n!  with  n  large)  to  the  familiar  Gaussian  form 
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P(8)  =  [2/irNf '[EXPHVsN)] 


(10) 


This  is  the  probability  that  a  particle  has  moved  s  mean  free  paths  in  N  collisions 
in  one-dimensional  motion.  Rewriting  this  in  terms  of  particle  speed  v  and  time  t,  and 
extending  to  three  dimensions,  yields 

P(r)  =  (2£/T;vt]^-^[EXP(-rV2vt£)]  (11) 


as  the  probability  of  diffusion  to  radius  r  from  the  origin,  in  time  t,  expressed  in 

one-dimensional  spherical  coordinates. 

#  - 

This  can  be  applied  in  two  ways  to  the  present  problem  of  diffusion  due  to 
numerous  small  angle  collisions.  One  way  is  to  determine  the  radial  displacement 
associated  with  a  given  number  of  collisions,  and  thus  to  obtain  a  measure  of  ion  dispersion 
and  associated  core  convergence  defocussing  due  to  same.  Another  is  to  estimate  the 
energy  change  due  to  small  angle  collisions  over  a  given  time  or  collision  number. 

This  latter  is  straightforward,  simply  noting  that  the  probability  P  equally  defines 
diffusion  in  energy  space,  where  r  ->  E  (E  is  total  energy  spread),  with  £  ->  6E^  as  the 
energy  change  per  collision  (the  "energy  mean  free  path"),  and  v  — >  f6E^  as  the  "speed"  in 
energy  space,  thus 

P(E)  =  (2/irft]^-®[EXP(-E^/2ft(6E^)^)]  (12) 

Here  the  energy  exchange  per  collision  is  just  6E^  =  Tr(Ze)  /Ljj,  where  Lj^  is  the 
Debye  length,  Lj^  =  (kT/4ime^)®'®,  as  usual,  or  =  (E^/6ime^)®’®  for  the  thermalized 
case  where  =  (3/2)kT. 
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The  mean  free  path  in  coordinate  space  is  just  £  =  l/nojj  where  <Tjj  is  the  Debye 
sphere  small  angle  collision  cross-section  (ff.  Lovberg<),  The  collision 

frequency  is  then  simply  f  =  v/£  =  (nvojj),  where  v  =  (2E^/m.)®'®.  Now  the  probability 
that  a  particle  will  scatter  into  some  energy  0  <  E  <  cd  must  be  exactly  unity.  Thus  the 
energy  probability  P(E)  can  be  normalized  to  C/P(E)dE  =  1,  integrating  from  0  to  cd,  from 
which  it  is  found  that  the  coefGcient  C  =  (irft/28E^),  so  that  the  probability  of  scattering 
to  an  energy  displacement  E,  from  initial  energy,  per  unit  energy,  is  just 

PJE)  =  lC2/irft)“'^l/6EJl[EXP(-EV2ft(6E_^)'‘)  (13) 


With  this  the  most  probable  value  of  the  scattered  energy  <E(t)>  at  any  time,  from 
small-angle  collisions  about  an  initial  energy  state  (i.e.  the  Gaussian  full-width  at 
half— maximum  amplitude)  is  found  from 


<E(t)>  =  [e  P^(E)  dE  =  (2ft/Tr)°’®(6E^) 

J  n 


(14) 


which  is  analogous  to  Lovberg’s  eq.  (12)  and  its  predecessor.  Using  the  Debye  collision 
cross-section  as  before,  this  gives  an  energy  dispersion  due  to  small— angle  scattering  over 
time  t  of 


<E(t)>  =  ir(2nvt)°-®(Ze)*  (15) 

For  10  keV  ions  uniformly  distributed  throughout  the  system  at  a  density  of 
n  =  lE12/cm  ,  and  with  a  system  lifetime  of  t  =  IE— 2  sec  (as  in  Lovberg’s  first  example), 
this  yields  <E(t)>  »  632  eV  or  ±  316  eV  as  the  deviation  from  mean  energy.  Note  that  the 
coUisional  time  assumed  gives  an  integrated  area  density  of  (ion)  collision  centers  along  the 
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A 

ion  path  of  about  lE18/cm  for  any  given  ion  transiting  the  uniform  system. 

As  previously  shown,  the  ion  system  transit  time  is  about  8tp^  —  Tr(R/v^),  where  v^ 
is  ion  speed  at  core  conditions  (assuming  a  relatively  flat  central  region  without  significant 
virtual  anode  formation).  Typically  the  transit  time  is  6t„  »  3.1E-6  sec  for  R  =  100  cm 
and  =  1E4  eV;  this  corresponds  to  about  G.  »  3.2E3  ion  transits  through  the  system. 
This  is  roughly  equivalent  to  an  electron  current  recirculation  factor  (outside  the 
stagnation  radius  r^  of  «  G.(m./m^)®‘®  »  1.4E5. 

The  uniform  distribution  of  ions  assumed  above  is,  of  course,  not  realistic  for  the 

A 

system.  Rather,  the  ion  distribution  will  vary  approximately  as  (1/r  )  outside  a  core  of 
radius  r^,  and  will  increase  faster  within  the  core.  With  this  distribution,  and  assuming  a 
3:1  density  jump  at  the  core  boundary,  as  used  for  convenience  in  earlier  calculations  (e.g. 
DTI/1989^,  Lovberg<),  the  integrated  value  of  <nvt>  «  2G.Jn(r)dr,  integrated  over  the 
range  0  <  r  <  R,  is  about  3.32(G.nj^).  For  the  G.  value  above,  this  gives  a  much  smaller 
value  for  collisional  density  than  before;  here  <nvt>  «  l.lE16/cm^  which  yields  an  energy 
spread  of  <E(t)>  =  63  eV  or  a  half-width  of  ±  31  eV. 

As  a  final  example  consider  a  core  density  of  n^  =  3E17/cm^,  a  core  ion  energy  of 
E^  =  3E4  eV,  r^  =  1  cm  and  the  same  number  of  ion  transits  as  in  the  previous  case. 
Then,  using  the  density  distributions  dted  by  Lovberg,  the  total  ion  path  density  per  unit 
area  is  <nvt>  =  3.92E21/cm*.  For  this  system,  the  energy  spread  becomes  (eq.  15) 
<E(t)>  «  2.82E4  eV  =  28.2  keV,  for  a  half-width  of  14.1  keV.  This  spread  seems 
acceptable  over  such  a  lifetime,  even  though  the  system  is  significantly  (but  not 
completely)  Maxwellianized. 

For  such  a  quasi-Maxwellianized  core-ion  thermal  energy  (temperature)  of 
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=  30  keV,  the  initial  radial  ion  injection  energy  must  be  about  (3/2)E^  or  45  keV.  The 
mean  upscatter  energy  dispersion  is,  as  above,  about  *  14.1  keV.  Allowing  for  2.5  times 
the  mean  ion  temperature  gives  a  well  depth  of  75  keV  which  is  sufGdent  to  trap  the  mean 
upscattered  ions.  Ions  inserted  into  the  well  at  the  45  keV  point  would  be  well  inside  the 
rippled  B  field  region,  and  less  prone  to  transverse  momentum  generation  than  for  higher 
energy  injection  at  larger  radii.  Lesser  well  depth  could  also  be  used,  but  at  a  price  in  ion 
loss  rate,  as  previously  described. 

It  is  important  to  note,  in  this  cormection,  that  the  upscatter  energy  dispersion 
given  by  eqs.  (14,15)  and  used  above  is  only  a  mean  value,  and  that  the  actual  upscatter 
energy  follows  the  Gaussian  distribution  of  eq.  (12).  This  is  another  way  of  representing 
the  eventual  stable  Maxwellian  distribution  that  will  obtain  for  long-lived  trapping  of  core 
ions. 


Note,  also,  that  the  fusion  power  output  of  the  last  example  system  will  be 
n^<CTV>E^(4Tt/3)(r^)^F(<r^>),  where  F(<r^>)  =  1  +  (l-3<r^>)/3  k  1.32.  For 
E^  =  30  keV,  <av>  «  5.7E— 16  cm  /sec  (Maxwellian  average),  the  fusion  power  generation 
becomes  «  8.0E8  watts  or  800  Mw.  Thus,  operation  with  DT  at  conditions  useful  for 
fusion  power  production  seem  tractable  within  the  constraints  of  collisional  upscattering. 

Of  course,  well  depths  must  be  chosen  large  enough  to  avoid  bulk  two—  and  three— body 
collisional  problems  in  all  events,  as  these  are  both  strongly  well-depth  dependent. 

As  a  last  topic,  consider  the  coordinate  spreading  introduced  by  small  angle 
scattering.  The  coordinate  probability  given  by  eq.  (11)  can  be  normalized  to  unity,  as  was 
done  for  spreading  in  energy  space,  by  integration  over  0<  r  <».  Noting  that  the  mean 
&ee.path  is  £  =  carrying  out  the  integration  gives  the  normalizing  coefficient  as 

K  s  (‘irvt)(n7Q/2)^,  and  the  probability  of  reaching  a  displacement  r,  per  unit  path  length. 
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becomes 


P.(0  =  [i>(ig/2iPrt]“‘(EXP(-(n<rj,/2vt)r2)l  (16) 

The  most  probable  displacement  over  any  given  time  period  <8r(t)>  is  just  the 
function  rP^(r)  integrated  over  the  complete  range  of  r.  This  yields 

<6r(t)>  =  (vt/2'Knajj)°'^  =  (6e^/ir)°®/(2mE^)®-^®  (17'i 

for  and  =  (E^/dime^)^'®  as  before.  Using  the  previous  examples,  the  mean 

displacement  due  to  small  angle  coUisional  distortion  of  the  initial  energy  distribution  is 
found  to  be  only  <6r(t)>  =  3.8E-3  cm  for  n^  =  lE12/cm^,  E^  =  1E4  eV,  and  t  =  lE-2 
sec.  For  the  more  realistic  case  of  non-uniform  density,  where  vt  =  2G.R  »  0.62E6  cm 
(rather  than  0.62E8  cm  as  in  the  initial  case)  this  is  found  to  be  only  8.6£^  cm  for  the 
same  ion  recirculation  factor  G..  It  thus  appears  that  energy  exchange  collisions  will  not 
contribute  directly  in  any  significant  way  to  core  convergence  spatial  broadening. 

However,  the  upscattering  of  energy  will  cause  ions  to  orbit  to  larger  radii,  and  thus 
to  enter  regions  of  larger  B  field,  in  which  the  B  field  structure  is  more  pronounced  (as  cf. 
to  near-sphericity  in  the  core  region).  This  (outward)  penetration  of  stronger  and  more 
rippled  B  fields  will  result  in  the  transformation  of  a  larger  fraction  of  ion  radial  energy 
into  transverse  momentum  than  for  orbits  to  smaller  radii.  This  larger  transverse 
momentum,  in  turn,  can  defocus  the  ions  in  their  subsequent  transits  of  the  core  region, 
and  may  even  result  in  the  trapping  of  ions  in  orbit  paths  well  outside  the  original  core 
radius. 

«  « 

The  upscattering  process  thus  results  in  core  convergence  broadening  from  the 
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indirect  effect  of  increased  transverse  momentum  generation  due  to  ion  orbit  expansion.  In 
the  absence  of  restoring  focussing  effects,  this  process  can  continue  only  until  the  ion  core 
offset  radius  is  at  a  distance  approximately  equal  to  half  of  the  gyro  radius  of  an  ion 
calculated  as  with  core  energy  in  the  cusp  surface  field,  at  which  time  the  ion  will  cease  to 
occupy  a  repetitive  orbit,  but  will  leave  the  system.  Such  restoring  effects  may  be  found  in 
transverse  collisions  in  the  core*  and/or  at  the  outer  radius  of  ion  motion,  as  described  by 
Rosenberg  and  Krall*  Further  exploration  of  this  coupling  of  energy  upscattering  to  ion 
orbit  distortion  and  momentum  transformation  is  required  to  obtain  a  good  understanding 
of  these  phenomena,  and  of  their  effects  in  various  systems. 
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